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Abstract 


Simultaneous wireless information and power transfer (SWIPT) provides a promising 
solution for enabling perpetual wireless networks. As energy efficiency (EE) is an im¬ 
portant evaluation of system performance, this thesis studies energy-efficient resource 
allocation algorithm designs in SWIPT systems. We first investigate the trade-off between 
the EE for information transmission, the EE for power transfer, and the total transmit 
power in a basic SWIPT system with separated receivers. A multi-objective optimization 
problem is formulated under the constraint of maximum transmit power. We propose 
an algorithm which achieves flexible resource allocation for energy efficiencies maxi¬ 
mization and transmit power minimization. The trade-off region of the system design 
objectives is shown in simulation results. Further, we consider secure communication in 
a SWIPT system with power splitting receivers. Artificial noise is injected to the com¬ 
munication channel to combat the eavesdropping capability of potential eavesdroppers. 
A power-efficient resource allocation algorithm is developed when multiple legitimate 
information receivers and multi-antenna potential eavesdroppers co-exist in the system. 
Simulation results demonstrate a significant performance gain by the proposed optimal 
algorithm compared to suboptimal baseline schemes. 
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Glossary 


Abbreviations 

AWGN 

Additive White Gaussian Noise 

GSIT 

Ghannel State Information at the Transmitter 

EE 

Energy Efficiency 

EH 

Energy Harvesting 

EH-EE 

Energy Harvesting Efficiency 

IR-EE 

Information Receiving Energy Efficiency 

MISO 

Multiple-Input Single-Output 

MIMO 

Multiple-Input Multiple-Output 

MOO 

Multi-objective Optimization 

MOOP 

Multi-objective Optimization Problem 

QoS 

Quality of Service 

RE 

Radio Erequency 

SDP 

Semi-definite Programming 

SIG 

Successive Interference Gancellation 

SINR 

Signal-to-Interference-plus-Noise Ratio 

SWIPT 

Simultaneous Wireless Information and Power Transfer 

TDD 

Time Division Duplexing 

WPT 

Wireless Power Transfer 


Operators 

(•)^ Hermitian transpose 

(•)* Complex conjugate 

0 All-zero matrix 

£{•} Statistical expectation 

Tr(-) Trace of a square matrix 

|j Absolute value 

II j| Euclidean norm 

det(-) Matrix determinant 

Null(-) Orthonormal null space of a matrix 

[x]"*" max{0,x} 

C^xM -phe space of all x M matrices with complex entries 

^^(m, S) A complex Gaussian random variable vector with mean vector m and 

covariance matrix S 



Glossary 


Symbols 

Bandwidth [Hz] 

Ai-j. Number of antenna equipped at the transmitter 

Nr Number of antenna equipped at the receiver 

Antenna power consumption [J/s] 

Static circuit power consumption [J/s] 

Total power consumption [J/s] 

^ Power amplifier efficiency 

r\ Energy conversion efficiency 
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Chapter 1. 
Introduction 


Wireless communication networks has been rapidly developing over decades. High 
speed and ubiquitous service is consistently expected in the evolution. This leads to 
tremendous energy demand for supporting the system operation. However, mobile 
devices with limited battery supply creates the bottleneck in providing continuous 
communication services. In particular, the slow improvement of battery capacity has 
hindered the fulfillment of high quality of service (QoS) requirements. Consequently, 
energy harvesting (EH) provides a new paradigm that enables self-sustainability for 
energy constrained wireless devices. Among different EH technologies, a promising one 
is wireless power transfer ^ (WPT) where communication terminals harvest energy from 
radio frequency (RE) signals. Recently, simultaneous wireless information and power 
transfer (SWIFT) has drawn much attention as an interesting and challenging scenario, 
since RE signal is a carrier of information and power concurrently. 

On the other hand, communication security is a critical issue as wireless commu¬ 
nications has become an indispensable media by which people may exchange secret 
information. As an alternative to traditional cryptographic techniques, physical layer 
security derives perfectly secure communication by exploiting the physical properties of 
wireless channel. 

In this chapter, we first give a brief overview of SWIFT and physical layer security. 
Then, we state the motivation of this thesis. 

1.1. Simultaneous Wireless Information and Power 
Transfer 

The integration of EH capability into wireless devices is a promising solution for prolong¬ 
ing the lifetime of communication networking. In practice, natural energy sources such 

^In this thesis, normalized unit energy is considered, i.e., Joule-per-second, which means the terms 
“energy" and “power" are interchangeable here. 
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as solar, wind, and geothermal are exploited. However, the challenge is that these re¬ 
newable sources are usually weather and location dependent which may not be suitable 
for portable mobile terminals. Fortunately, WPT is a promising solution since it allows 
energy harvester scavenge energy from relative stable and controllable electromagnetic 
waves in both indoor and outdoor environment. In particular, RF signal as a carrier 
of energy is an abundant source for WPT. Nowadays, EH circuits are able to harvest 
microwatt to milliwatt of power over the range of several meters for a transmit power 
of 1 Watt and a carrier frequency less than 1 GHz [1]. Thus, RF energy can be a viable 
energy source for devices with low-power consumption, e.g. wireless sensors [2, 3]. 
Moreover, RF EH provides the possibility for simultaneous wireless information and 
power transfer (SWIPT) [4, 5]. As a carrier of both information and energy, RF signal 
has a remarkable superiority that unifies information transmission and power transfer, 
and embrace energy saving by enabling tremendous energy consumed by wireless signals 
recyclable. As a result, SWIPT has gained recent attention in academic research. 

In a SWIPT system, when information signals conveyed from the transmitter to the 
intended receivers via wireless channels, energy harvesters (the same or other receivers) 
can harvest energy from the information signals due to the broadcasting nature of 
wireless channels. Different from conventional wireless communication systems where 
data rate is the most fundamental system performance metric, the amount of harvested 
energy has become an equally important QoS requirement in SWIPT systems. Thus, 
new resource allocation algorithms are needed to fulfill the emerging need [4]-[20]. 
In [4]-[7], the fundamental trade-off between channel capacity and harvested energy 
was studied for frequency flat fading channel and frequency selective fading channel. 
Specifically, an ideal receiver managing synchronous information decoding and energy 
harvesting from the same received signal was assumed in [4] and [5]. However, the 
signal used for information decoding, cannot be reused for EH due to the limitation 
of current practical circuits. Subsequently, the author in [6] and [7] proposed three 
different t5q5es of receivers, namely, power splitting, separated, and time switching 
receivers. In particular, a power splitting receiver, cf. Eigure 1.1, splits the received 
signal into two power streams with a power splitting ratio p, i.e., p portion of the 
received signal remains for information decoding, and the other 1 — p portion is reserved 
for EH. The power splitting scheme also generalizes the separated receivers scheme, i.e., 
independent information decoding is performed when p = 1 and independent EH is 
done when p = 0. A time switching receiver means switching time slots for information 
decoding and EH successively. 

Furthermore, [8] and [9] focused on transmit beamforming design in multiple-input 
single-output (MISO) SWIPT systems for separated and power splitting receivers, re- 
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Figure 1.1.: Power splitting receiver. 


spectively. In [10], beamformers were optimized for the maximization of sum harvested 
energy under the minimum required signal-to-interference-plus-noise ratio (SINK) con¬ 
straints for multiple information receivers. In [11], [12], and [13], energy-efficient 
SWIPT, was studied in multi-carrier systems, where power allocation, user and subcarrier 
scheduling were considered. It was shown that EE could be improved by implement¬ 
ing SWIPT. In [14], a transmission strategy was proposed for multiuser multiple-input 
multiple-output (MIMO) SWIPT systems with interference channel. Recently, multiuser 
scheduling, which exploits multiuser diversity for improving the system performance 
of multiuser systems was considered for SWIPT in [15] and [17]. On the other hand, 
SWIPT also envisions new opportunities for cooperative communications. In [19], the 
performance of SWIPT systems was analyzed for different relaying protocols. For cogni¬ 
tive radio networks, [20] focused on the cooperation between primary and secondary 
systems at both the information and energy levels. In these literatures, SWIPT demon¬ 
strates significant gains in many aspects, for instance, energy consumption, spectral 
efficiency, and time delay. 


1.2. Physical Layer Security 

Nowadays, wireless communication security has become an extremely important issue. 
Traditionally, communication security relies on crjqjtographic technologies that are 
applied in the application layer of wireless networks. These algorithms usually require 
reliable key distribution and may encounter high computational complexity. As an 
alternative, physical layer security is a viable solution for ensuring communication 
security. It guarantees secure communication by exploiting the physical properties of 
wireless channels. In [21], the basic idea of physical layer security was first proposed, 
cf. Figure 1.2. It is verified that confidential messages can be reliably exchanged 
between a transmitter (Alice) and a receiver (Bob) if the receiver enjoys a better channel 
than an eavesdropper (Eve). In particular, by exploiting the extra degrees of freedom 
offered by multiple transmitting antennas, the channel between the transmitter and the 
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Noise 



Figure 1.2.: A general wiretap model. 


eavesdropper can be weakened by injecting a properly designed artificial noise. Thus, 
perfect communication security is achievable. 

Under this context, secure communication guaranteed by artificial noise generation 
led to a series of discussion. In [22] , a power allocation algorithm was proposed for 
the maximization of ergodic secrecy capacity. In [23], the author maximized secrecy 
capacity in a system with multiple single-antenna eavesdroppers. [24] studied secure 
OFDMA systems by energy-efficient resource allocation algorithm design. A trade-off 
between energy efficiency and secure communication was revealed. Further, [25] 
focused on secure resource allocation and scheduling in OFDMA relay networks, where a 
passive multi-antenna eavesdropper was considered. In [26], the author studied secure 
communication in cognitive radio networks. A robust resource allocation algorithm was 
proposed to ensure video communication secrecy in the secondary system. 

In fact, a large amount of power is allocated to artificial noise for providing secure 
communication. This implies that artificial noise can act as a potential energy source 
for WPT while ensuring secure communication at the same time. Resource allocation 
algorithm design for secure communication with SWIFT was studied in [27]-[35]. In 
[27], an additive energy signal was adopted to facilitate power transfer and to ensure 
secure communication for separated information and energy receivers. [28], [30], and 
[32] studied power-efficient resource allocation algorithm for secure SWIFT systems 
by minimizing the total transmit power. In [28], robust beamforming was studied 
in a system with imperfect CSIT, where both artificial noise and energy beams are 
used in order to provide secure communication and to improve WPT. Besides, in [29] 
multi-objective optimization (MOO) approach was applied to jointly optimize multiple 
system design objectives in a secure communication system with SWIPT. Especially, 
EE of energy harvesting was maximized. In [35], SWIPT in secure communication 
system was combined with cognitive radio, where multiple objectives, including total 
transmit power minimization and EH efficiency maximization were considered in a 
multi-objective optimization problem (MOOP). 
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1.3. Motivation 

As aforementioned, in SWIPT systems, the system performance is evaluated on both 
information delivery and power transfer. A trade-off naturally arises when considering 
resource allocation on both aspects. The trade-off between data rate and the amount 
of harvested energy has been investigated in recent literatures. However, the trade-off 
between EE of information delivery and EH has not been considered so far. On the 
other hand, the conflicting system design objectives, i.e., information receiving EE maxi¬ 
mization, energy harvesting efficiency maximization, and transmit power minimization, 
leads to multiple resource allocation algorithm designs. To provide a resource allocation 
algorithm that can flexibly achieve multiple objectives, multi-objective system design for 
SWIPT requires to be studied. 

Furthermore, communication security is a serious issue in SWIPT systems. In particular, 
when the transmit power of the information signal is increased to facilitate SWIPT, 
the signal becomes more vulnerable to eavesdropping due to a higher potential for 
information leakage. Thus, communication security arises as a new QoS concern in 
SWIPT systems. On the other hand, artificial noise can serve as a energy source for 
WPT and be harvested to extend the lifetime of power-constrained devices. In terms 
of power-efficient resource allocation in secure SWIPT systems, recent literatures limit 
the system configuration with a single information receiver and multiple single-antenna 
eavesdroppers. However, optimal resource allocation for secure communication in SWIPT 
systems with multiple desired power splitting receivers and multi-antenna potential 
eavesdroppers remains an unsolved problem. 

The rest of the thesis is organized as follows. In Chapter 2, we study multi-objective 
SWIPT design for separated receivers. The energy efficiencies for information trans¬ 
mission and EH are investigated by taking into account the maximum transmit power 
constraint. A resource allocation algorithm considering the trade-off between multiple 
objectives is proposed. In Chapter 3, we focus on SWIPT in secure communication 
systems with multiple power splitting receivers and multi-antenna eavesdroppers. A 
power-efficient resource allocation algorithm is proposed. Finally, we summarize the 
contributions of this thesis in Chapter 4. Appendix A contains basic theories of optimiza¬ 
tion problem. Appendix B includes the proofs of the theorems and propositions in the 
thesis. 
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Chapter 2. 

Multi-Objective SWIPT with 
Separated Receivers 


In this chapter, we study multi-objective power allocation for energy-efficient SWIPT 
with separated receivers. In a MISO system, information signals and energy beams 
are transmitted simultaneously to jointly support information delivery to a information 
receiver and energy supply to a energy harvester. Under a maximum transmit power 
constraint, we focus on three desired system design objectives, namely, information 
receiving efficiency (IR-EE) maximization, energy harvesting efficiency (EH-EE) maxi¬ 
mization, and total transmit power minimization. In particular, we jointly optimize the 
information beamforming vector and covariance matrix of the energy signal to achieve 
the considered system objective. The problem is formulated as a non-convex MOOR To 
deal with the fractional objective functions, Charnes-Cooper transformation method is 
adopted. Subsequently, the transformed problem is solved by semi-definite program 
(SDP) relaxation approach. We prove that the SDP relaxation is tight. In particular, 
a tractable structure of the optimal solution is verified. Simulation results shows the 
trade-off between IR-EE, EH-EE, and the total transmit power. 


2.1. System Model 

We focus on a downlink MISO system with SWIPT. The system consists of one multi¬ 
antenna transmitter, one single-antenna information receiver, and one single-antenna 
energy harvester. The transmitter is equipped with antennas. It sends precoded infor¬ 
mation signal and energy beams simultaneously to facilitate information transmission 
and power transfer, cf. Eigure 2.1. The transmission is divided into time slots. The 
transmitted signal in each time slot is given by 

x = WjS -I- Wg, 


( 2 . 1 . 1 ) 
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Transmitter 





Energy harvester 

Figure 2.1.: A SWIPT system with separated receivers. 


where s e C is the information-bearing symbol with £||s l^j- = 1. Wj e c^txi is the 
corresponding precoded beamforming vector for the information receiver, w^ e is 

the energy signal beamforming vector facilitating energy transfer to the energy harvester. 
The energy beamforming vector w^ is modeled as a complex Gaussian pseudo-random 
sequence as w^ ~ W^), where = E^w^Wg | is the covariance matrix of the 

energy signal. Assume w^ is generated at the transmitter by a pseudo-random sequence 
generator with a predefined seed. The seed can be delivered to the information receiver 
before effective information transmission. Thus, the interference of energy signal can be 
totally cancelled at the information receiver. 

We assume a narrow-band slow fading channel between the transmitter and receivers. 
The channel is assumed to be perfectly known at the transmitter. Then, the received 
signals at the information receiver and energy harvester are expressed as 

yiR = h"(wiS-FwE)-Fni, (2.1.2) 

Jeh = g"(wiS + We) + Ue, (2.1.3) 

where h e channel vector between the transmitter and the information 

receiver, and g e channel vector between the transmitter and the energy 

harvester. They capture the joint effect of multipath fading and path loss. U; e C and 
Ue e C are additive white Gaussian noise (AWGN) at the information receiver and energy 
harvester, respectively, which are distributed as ‘<^^(0, cTj^) and ‘i^^(0, cr^). 

Information receiver focuses on decoding the information signal. The achievable rate 
(bit/s/Hz) at the information receiver can be described as 

R = log2 (l + = log2 (^1 + ^wfHwi), 


(2.1.4) 
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where H = hh^. At the same time, both the information signal and the energy signal 
can act as RF energy source for the energy harvester due to the broadcast nature of 
wireless channels. According to the law of energy conservation, the harvested energy 
is proportional to the received signal power. The total harvested energy at the energy 
harvester is given by 

f^harv(wi,WE) = i7(^|g"wiP + |g"wEp) = Tj (^wf Gwj + Tr(GWE)), (2.1.5) 

with G = gg^. 17 is the energy conversion efficiency, which is a constant with 0 < 17 < 1. 
It implies a energy loss in the process of converting the received RF energy to electrical 
energy for storage. We ignore the thermal noise at the receiving antenna as it is relative 
small compared to the received signal power. 


Apart from system throughput, EE is also a fundamental system performance metric 
in modern communication networks. EE is generally defined as the ratio between 
system throughput and total power consumption. We first model the total power 
consumption (Joule-per-second) by taking into account the transmit power consumption 
and additional hardware power dissipation at the transmitter which can be described as 

||wt|P + Tr(WE) 

We) = -^- +iVTPant + Pc- (2-1.6) 

^ is the power amplifier efficiency, which is a constant with 0 < <^ < 1. The first term in 
( 2 . 1 . 6 ) is the total power consumption in the power amplifier. accounts for the 

dynamic circuit power consumption proportional to the number of transmitting antenna. 

denotes the power dissipation at the transmitting antenna, including the transmit 
filter, mixer, frequency S 5 mthesizer, digital-to-analog converter (DAG), etc. denotes 
the fixed circuit power consumption due to baseband signal processing. 


Based on the general concept of efficiency, we define IR-EE and EH-EE as 


^ir(wi,We) = 

R 

log 2 (l + ;^wfHwi) 

_ 

(2.1.7) 

P 

^tot 

“ (||Wi|P + Tr(WE))/?+iVTPant + J’c 

$eh(wi,We) = 

P 

■^harv 

^(w^Gwi + Tr(GWE)) 

( 2 . 1 . 8 ) 

p 

^tot 

“ (||Wi||2 + Tr(WE))/?+iVTPant + f^c’ 


respectively. 
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2.2. Problem Formulation 

In SWIPT system, IR-EE maximization, EH-EE maximization, and total transmit power 
minimization are all desirable for system design. In this section, we first propose three 
problem formulations for single-objective system design for SWIPT. Each single-objective 
problem describes one important aspect of the system design. Then, we consider the 
three system design objectives jointly by MOO. 

The first system design objective is the maximization of IR-EE. The optimization 
problem is formulated as 

Problem 2.1. IR-EE Maximization: 

maximize $jr(wi, W^) 

WeSH^t^w, 

subject to Cl : Hwjp + TrCWg) < 

C2: We^O. 

The second system design objective is the maximization of EH-EE. The problem 
formulation is given as 

Problem 2.2. EH-EE Maximization: 

maximize 

WeSH^t^w, 

subject to Cl, C2. 

The third system design objective is the minimization of the total transmit power at 
the transmitter. The problem formulation is proposed as 

Problem 2.3. Total Transmit Power Minimization: 

minimize 11 Wj | ^ -I- Tr (W^) 

WEeH~T,w, 

subject to Cl, C2. 

As the above problem formulations stated, IR-EE, EH-EE, and the total transmit power 
are independently optimized respectively. In each single-objective problem, information 
beamforming vector, Wj, and the covariance matrix of the energy signal, W^, are jointly 
designed by considering the maximum transmit power constraint Cl. In addition, 
covariance matrix Wg should be a positive semi-definite Hermitian matrix as indicated 
in constraint C2. 

For the sake of notational simplicity, we denote the objective functions in the above 
problems as Fj(wi,We), j = 1,2,3. We note that Problem 2.3 is a trivial problem 
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with optimal value zero since the transmitter does not need to provide any QoS to the 
receiver. Yet, Problem 2.3 plays an important role in the following when we study the 
multi-objective power allocation algorithm design. Without loss of generality, Problem 
2.3 can be rewritten as an equivalent maximization problem in order to represent 
the three problems consistently. The corresponding objective function is written as 
F3(wi,WE) = -(||wJ|2 + Tr(WE)). 

In practice, these three independent optimization objectives are all desirable from 
the system operator perspective. However, there are non-trivial trade-off between 
them. In order to optimize these conflicting system design objectives systematically and 
simultaneously, we apply MOO, cf. Appendix A.l. 

A common approach to formulate a MOOP is via "a prior method". This method allows 
the system designer to indicate the relative importance of the system design objectives 
before running the optimization algorithm. In particular, a sequence of scalars, which 
is known as "preference parameters" or "weights", are a prior specified to scalarize 
system designer’s preference on different objectives. There are many scalarization 
methods. Here we adopt weighted min-max method [36]. As introduced in Appendix 
A.l, the optimal point(s) of a MOOP is defined by Pareto optimality. All Pareto optimal 
points, which form the Pareto optimal set, are important to the system designer. In 
fact, weighted min-max method can provide the complete Pareto optimal set by varying 
the preference parameters, despite the non-convexity of the MOOR Based on weighted 
min-max method, we incorporate three system design objectives into a MOOI] which is 
formulated as 


Problem 2.4. Multi-Objective Optimization Problem: 


minimize 

We£H~t,w, 

subject to 


~ F/wj, We))} 

Cl, C2, 


where F* is the optimal objective value with respect to Problem j. coj is the weight 
imposed on objective function j subject to 0 < <Uj < 1 and coj = 1, which indicates 
the system designer’s preference on jth objective function over the others. In extreme 
case, when coj = 1 and to j = 0, Vi 7 ^ j, Problem 2.4 is equivalent to the single-objective 
optimization problem j. 

In this MOOI] we investigate the complete trade-off region between the three objectives 
regarding to system power allocation. To this end, we only take maximum transmit 
power constraint into consideration. In case other QoS constraints are imposed into the 
MOOI] a smaller Pareto optimal set can be obtained, which is actually a subset of the 
complete trade-off region. 
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2.3. Multi-Objective Power Allocation Algorithm Design 


It can be observed that Problem 2.1 and Problem 2.2 are non-convex due to the fractional 
form of the objectives which leads to the non-convexity in Problem 2.4. In order to 
obtain a tractable solution, we first transform the non-convex objective functions by 
Charnes-Cooper transformation method. Then, the transformed problems are solved by 
SDP relaxation approach. 

We first reformulate aforementioned three single-objective optimization problems by 
defining a set of new optimization variables as follows: 

Wi = wiwf, 9 = \ Wi = 0Wi, and (2.3.1) 

Ptot(Wi,WE) 

Then the original problems can be rewritten with respect to the new optimization 
variables {Wj, Wg, 6], which are given by 


Problem 2.5. Transformed IR-EE Maximization Problem: 


m^imize 

w,,WEeH"T,e 

subject to 


01 og2(l + 


Tr(HWi) 

6(j^ 


) 


a : Tr(Wi + WE)<0P^ax, 

^ : Wj ^ 0 , We ^ 0 , 

^ : Rank(Wi) < 1, 

_ Tr(Wj + We) 

C4: -^- ^- + 6iN^P,,, + P,)<l, 


C5 : e>0. 


Problem 2.6. Transformed EH-EE Maximization Problem: 

maximize ri Tr(G(W, -I- Wp)) 

W,,WE£H~T,e 

subject to Cl — C5. 


Problem 2.7. Transformed Total Transmit Power Minimization Problem: 


m^imize 

WpWEeH~T,0 

subject to 


-^i^-N^P,n,-Pc) 
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Denote the transformed objective function as Fj(Wi,We, 6), j = 1,2,3. Constraints 
Wj ^ 0, Wj e and Rank(Wi) < 1 are imposed to guarantee that Wj = 0WjW^. 
Constraints C4 and C5 are introduced due to the proposed transformation. 

Furthermore, in order to simplify the following algorithm design, we first normalize 
the transformed objective functions due to their different ranges and dimensions. A 
robust transformation, regardless of the original range or dimension of the objective 
function, is given as follows [36], 


-F/Wi,We,0)-f/ 

Fj (Wi, We, 6) = _o - 

Fj -Fj 


(2.3.2) 


—^ —0 

where Fj and Fj are the maximum and minimum value of the jth transformed objective 
function, i.e., Fj° < Fj(Wi,We, 0 ) < Fj . Fj can result from the transformed single¬ 
objective problems. Then, the transformed objective functions are normalized to range 
[ 0 , 1 ]. 

Regarding to the MOOP 2.4, the objective function can be rewritten in a normalization 
form as max^^^ 2,3 ~ 0))}. A common approach for handling such a 

min-max optimization problem is to introduce an auxiliary optimization variable. Then, 
the MOOP can be transformed into its equivalent epigraph representation [37], which is 
given by 


Problem 2.8. Transformed MOOP: 


jmnimize t 

W,,WE£H“T,e,T 

subject to Cl — C5, 

^ : cojil - f:"“'(Wi, We, 0)) < t, Vj, 

where t is the auxiliary optimization variable. We note that the optimal value of Problem 
2.8 lies between zero and one. 

Now, we introduce the following proposition. 

Proposition 2.1. The transformed problems 2.5-2.8 are equivalent transformations of the 
original problems 2.1-2.4, respectively. 

Proof. Please refer to Appendix B.l. □ 

Based on Proposition 2.1, we can recover the solution of the original problems based 

—^ —0 

on (2.3.1). In particular, the optimal value Fj or the lower bound Fj of the jth 
transformed objective function equal to that of the jth original objective function, i.e.. 
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Fj = Fj, Fj^ = F°, j = 1,2, 3. Thus, in (2.3.2), are trivial results from the original 

Q Q Q ^ 

objective functions, that are F^ = f 2 = Oj ^3 = ^max- ^l^o denote Fj simply as 

j;* = 4>*^, j;* = and r,* = o. 

We note that if Problem 2.8 can be solved optimally by an algorithm, then the 
algorithm can also be used to solve Problem 2.5-2.7, since Problem 2.8 is a generalization 
of Problem 2.5-2.7. Thus, we focus on the method in solving Problem 2.8. It is evident 
that Problem 2.8 is non-convex due to the rank-one beamforming matrix constraint 
C3 : Rank(Wi) < 1. Now, we apply the SDP relaxation by removing constraint C3 from 
Problem 2.8. As a result, the SDP relaxed problem is given by 

Problem 2.9. SDP Relaxed Transformed MOOP: 

jiunimize t 

Wi,WEeH"T.e,T 

subject to Cl, C2, C4, C5, C6, 

which is a convex SDP problem and can be solved by numerical convex program solvers 
such as CVX [38]. In particular, if the obtained solution Wj of the SDP relaxed problem 
satisfies constraint C3, i.e., Rank(Wj) < 1, then it turns out to be the optimal solution. 
Then, the optimal beamforming vector w* of the original problem can be achieved 
by solving the relaxed problem and recovering from the invertible mapping equations 
(2.3.1). Now, we study the tightness of the SDP relaxation by the following theorem. 

Theorem 1. The optimal solution of Problem 2.9 satisfies Rank(Wj) = 1 and Rank(Wj.) < 
1. In particular, an optimal solution with Rank(Wj) = 1 and = 0 can always be 
constructed. 

Proof. Please refer to Appendix B.2. □ 

Therefore, the adopted SDP relaxation is tight. Besides, Problem 2.5-2.7 can be solved 
by SDP relaxation as solving Problem 2.9. 

-^ -* 

Next, we construct an optimal solution with Rank(Wj) = 1 and Wg = 0 based on 
Theorem 1. We redefine the optimization variable W; as 

Wj = Auu^, u = [ui, U 2 ,..., UjvJ and Wg = 0, (2.3.3) 

where u e ^ orthonormal vector, i.e., |“i|^ = 1- According to (2.3.1) 

and (2.3.3), we have Wj = Then, MOOP 2.8 can be reformed with respect to the 
optimization variables {u. A, 0, t} as follows: 
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Problem 2.10. 

minimize 

u,A, 6 ,t 

subject to 


JVt 


Cl: xY,\ui\"<ep, 


i=l 


C2 : A > 0, C3 : 0 > 0, 




C4: 




+ 0(iVTPant + 0c)<l, 


— z 0 r 

C5: mifl-—log2M + 

^IR 


A 




C6 : <u, I 1 




** 

EH 


A 


Nt 




i=l 


0cr^ 


)<T, 


-))<T, 


C7: m3-^(^-iVTPant-0c)<^, 

mciY 


where h; and g;, i e {1,.. .,Nj}, are the elements of channel vectors h and g, respectively. 
In order to investigate the structure of vector u, we analyze the Karush-Kuhn-Tucker 
(KKT) conditions of Problem 2.10 by introducing the Lagrangian function. The La- 
grangian function is given by 


+ 


+ 


if (u, A, 0, T, /X, V, /Cl, K 2 , /C3, C) 


(2.3.4) 


Nt 


T + m(A 2 - 0P„T J + V + 0(iVTPant + 0c) - 1) - 


! = 1 




K, 


K. 


0 

^IR 


A 




Qcr^ 


■)) -^] 


CO. 


:(1 
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$ 


A 


EH 


Nt 

2 


) - t] + /C3 

i=l 





? 1 


’P '6 

^ max 


(g - iVT0ant-0c)-^], 


where ju, v, /c^, K 2 , k^, C are dual variables associated with the corresponding constraints, 
respectively. Constraint C2 is captured in the solution when deriving KKT conditions 
in the following. Since Problem 2.10 satisfies Slater’s constraint qualification and is 
convex with respect to the optimization variables, strong duality holds. Then, based on 
KKT optimality conditions, the gradient of Lagrangian function with respect to Uj, the 
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element of u, vanishes, from which we can result in 


U; 


where a 


co^ahi-\-co2bgi, 


and b = 


rv * V 




i=l 


Similarly, consider KKT condition with respect to A, which is given by 

-I + 


x = e 




<J 


(n(2)(M+|-^ SS.u-g, ) 


(2.3.5) 


(2.3.6) 


As we can see, (2.3.5) and (2.3.6) imply the structure of beamforming vector W; by 
considering W; = In particular, (2.3.5) indicates the direction of the information 

signal. (2.3.6) shows that the power allocation for the information signal follows the 
policy of water-filling solution. For specific case, when IR-EE is considered and EH-EE is 
discarded, i.e., 0 and CO 2 = 0 , information beamforming vector Wj is aligning to 

the direction of channel vector h according to (2.3.5). Since |“i|^ = obtain 


Wi= Tp 



where p = 


ln( 2 )(p + v/<^) 



(2.3.7) 


On the other hand, when IR-EE is not taken into account and EH-EE is maximized, 
i.e., = 0 and CO 2 7 ^ 0 , the beamforming vector Wj directs to the energy harvester by 

following the direction of channel vector g as (2.3.5) indicates. Especially, Problem 2.10 
becomes a linear programming with respect to A. In extreme case, if transmit power 
minimization is not considered either, i.e., co^ = 0 , we solve a single-objective problem 
for EH-EE maximization. Then, the optimal solution is given as 


w, 




g 

1 lgll’ 


EH 


’l^'maxllglP 


+ •^T^’ant + Pc 


(2.3.8) 


Furthermore, when both IR-EE maximization and EH-EE maximization are active 
objectives, i.e., 7 ^ 0 and CO 2 7 ^ 0, Wj is designed as a dual use beamforming vector for 

simultaneous information delivery and power transfer. (2.3.5) shows that it incorporates 
the directions of both channel vectors h and g. 
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2.4. Results 

In this section, we present simulation results to demonstrate the system performance of 
multi-objective system design. The simulation parameters are summarized in Table 2.1. 
In particular, we adopt the TGn path loss model [39] . The multipath fading is modeled 
as Rician fading with Rician factor 3 dB. Assume the carrier center frequency as 470 MHz 
with bandwidth 200 kHz. At the transmitter, we set the djmamic power consumption 
Pant = 75 mWatt per antenna, static circuit power consumption Pc = 1 Watt, and the 
power amplifier efficiency ^ — 0.4. The maximum transmit power is given as = 1 
Watt. Two receivers, namely, information receiver and energy harvester, are uniformly 
located between the reference distance 1 meters and the maximum service distance 
10 meters. Each receiver is equipped with a single antenna with antenna gain 10 dBi. 
Assume the noise covariances at the information receiver and the energy harvester are 
the same, i.e., erf = = cr^. We set cr^ = —47 dBm which includes thermal noise at 

a temperature of 290 Kelvin and signal processing noise. The signal processing noise 
is caused by a 12 -bit uniform quantizer employed in the analog-to-digital converter at 
the analog front-end of each receiver. At the energy harvester, the energy conversion 
efficiency for converting RF energy to electrical energy is 17 = 0.8. In this setting, 
multiple channel realizations are simulated, where both pass loss and multipath fading 
effects are taken into account. 


Table 2.1.: Simulation Parameters 


Carrier center frequency 

470 MHz 

Bandwidth 

m = 200 kHz 

Single antenna power consumption 

P,,, = 75 mW 

Static circuit power consumption 

Pc = 1 w 

Power amplifier efficiency 

d 

II 

Antennas gain 

10 dBi 

Noise power 

cr^ = -47 dBm 

Rician factor 

3 dB 

Reference distance 

1 meters 

Maximum service distance 

10 meters 

Energy conversion efficiency 

00 

d 

II 


In the following, we show the trade-off region between multiple system objectives 
from two aspects. In one aspect, we examine the trade-off between the average IR-EE, 
EH-EE, and transmit power in terms of system EE, which is shown in Figures 2.2, 2.4, 
2.5, 2.6, 2.7, 2.9, and 2.11. On the other hand, from the aspect of system throughput, 
the trade-off between average achievable rate, harvested energy, and transmit power 
is illustrated in Figures 2.3, 2.8, 2.10, and 2.12. For comparison, we also propose 
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Figure 2.2.: System performance trade-off region between IR-EE, EH-EE, and transmit 
power. 



Figure 2.3.: System performance trade-off region between achievable rate, harvested 
energy, and transmit power. 


a baseline scheme, where MOOP of achievable rate maximization, harvested energy 
maximization, and transmit power minimization is solved. The system performance are 
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compared between the proposed EE algorithm and the baseline scheme in Eigure 2.7- 
Figure 2.12. 

Figures 2.2 and 2.3 give the 3-dimension trade-off regions of the system energy 
efficiency and system throughput, respectively, for 8 transmitting antennas. The 3- 
dimension trade-off regions are obtained by solving the MOOP 2.4 with different sets of 
weights on the system design objectives. Specifically, the points consisting of the regions 
are calculated out by uniformly varying the weight coj with a step size of 0.04 such that 
cOj = 1. It can be observed in Figure 2.2 that the trade-off region between IR-EE, 
EH-EE, and transmit power is formed by the points gradually spreading from the right 
bottom corner to the left top corner. In particular, both IR-EE and EH-EE grow rapidly 
for small transmit power. When the transmit power is high, EH-EE continues increasing, 
however, IR-EE declines. On the other hand. Figure 2.3 illustrates that high transmit 
power supports the increment of both the achievable rate and harvested energy. 


X 10 ® 



Figure 2.4.: Trade-off region between IR-EE and EH-EE. 

In addition, for a better illustration, we also provide different side-views of the 2- 
dimension trade-off region in Figures 2.4, 2.5, and 2.6 for revealing the trade-offs 
between different pairs of objective functions. Figure 2.4 shows the trade-off between 
IR-EE and EH-EE. Figure 2.5 shows the trade-off between IR-EE and transmit power. 
Figure 2.6 shows the trade-off between EH-EE and transmit power. It can be observed 
from these figures that IR-EE and EH-EE are partially aligned with each other for small 
transmit power. In particular, IR-EE and EH-EE both increase rapidly when transmit 
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Figure 2.5.: Trade-off region between IR-EE and transmit power. 



Eigure 2.6.: Trade-off region between EH-EE and transmit power. 


power grows from zero. However, IR-EE reduces dramatically in the high transmit 
power regime which results in bell-shaped curves as shown in Figures 2.4 and 2.5. This 
diminishing return of IR-EE is due to the slow logarithmical growth of the achievable 
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rate in the high transmit power regime while the transmit power is linearly increasing. 
In contrast, EH-EE is monotonically increasing with the increasing transmit power as 
shown in Eigure 2.6. In other words, more energy is carried by the transmit signal, 
higher EH-EE can be achieved. This is thanks to the linear relationship between the 
harvested energy and the transmit power. 

Notably, the trade-off region in Figure 2.5 is non-convex. In other words, the proposed 
multi-objective system design algorithm is able to obtain the non-convex feasible region, 
despite of the non-convexity of the MOOR Besides, the three extreme points in Figure 2.4 
correspond to the three single-objective functions, respectively. The zero point for IR-EE 
and EH-EE in Figure 2.4 corresponds to the zero transmit power in Figures 2.5 and 2.6, 
which represents the minimum transmit power. It is the optimal value of single-objective 
Problem 2.3 which can also be obtained by solving the MOOR with CO 3 = 1. The second 
extreme point in the middle of the curve in Figure 2.4 is the maximum IR-EE, i.e., the 
optimal value of single-objective Problem 2.1, which can also be obtained by solving the 
MOOR with = 1. The third extreme point at the tail in Figure 2.4 demonstrates the 
maximum EH-EE, which is the optimal value of single-objective Problem 2.2. It can also 
result from the MOOR with (O 2 = 1. 

In Figures 2.7 and 2.8, the average IR-EE versus the average EH-EE, and the average 
achievable rate versus the average harvested energy are showed, respectively. These 
curves are obtained by solving the MOOR for co^ — 0 and 0 < coj < l,j = 1,2, where 
the value of cOj is uniformly varied with a step size of 0.01 such that coj = 1. Figure 
2.7 shows the trade-off between IR-EE and EH-EE when the objective of transmit power 
minimization is not considered. We can see that IR-EE is monotonically decreasing as 
EH-EE increasing, since the objective preference shifts from IR-EE to EH-EE, i.e., co^ 
decreases and CO 2 increases. Interestingly, we have a distinct dropping point at the tail of 
the curve corresponding to = 0 and C 02 = 1. This point indicates the solution of the 
single-objective problem of EH-EE maximization. Based on Theorem 1 and Appendix B.2, 
we have Rank(WE) = 1 at this point instead of = 0 in other points. In other words, 
the energy signal occupies a part of the total available power. Thus, IR-EE drops due to 
a smaller power allocation on information signal. Moreover, compared to the baseline 
scheme, it is obvious that IR-EE of the proposed EE algorithm achieves a significant gain. 
Besides, when the number of transmitting antenna is increased from Nj = 4 to N-^ = 8, 
the trade-off region is enlarged in both EE algorithm and the baseline scheme. Since 
extra degrees of freedom offered by more transmitting antennas can be exploited. Thus, 
the system performance on EE is improved. However, IR-EE for Aiy = 8 is smaller than 
that for = 4 in the low transmit power regime. This can be explained that with small 
transmit power, the achievable rates for both Nj = 8 and Nj = 4 are quite small. But the 
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total power consumption is large for Nj = 8 due to a high antenna power dissipation, 
which result in a smaller IR-EE. 


x10® x10® 




Eigure 2.7.: Average IR-EE versus 
average EH-EE. 


Eigure 2.8.: Average achievable rate versus 
average harvested energy. 


In terms of of system throughput, we see from Figure 2.8 that there is a consistent 
trend between the achievable rate and the harvested energy for EE algorithm. This 
is because as EH-EE is gradually emphasized, the amount of harvested energy grows 
due an increment of transmit power. On the other hand, according to Theorem 1 and 
Appendix B.2, Rank(Wj) = 1 and = 0 when IR-EE and EH-EE are both optimized. It 
means that the information signal occupies all the available transmit power. Thus, the 
information signal becomes stronger with the increasing transmit power which brings 
the improvement of achievable rate. As a result, the alignment between the achievable 
rate and the harvested energy occurs. In contrast to this consistency in the proposed 
algorithm, the achievable rate and the harvested energy in the baseline scheme conflict 
with each other. In particular, high data rate corresponds to low harvested energy, and 
vice versa. Moreover, it is noted that the curve stretches to the right end with a distinct 
drooping point as in Figure 2.7. This is caused by the same reason aforementioned. 
Besides, when more transmitting antennas are equipped at the transmitter, the system 
throughput with respect to the achievable rate and the harvested energy are both 
improved since extra degrees of freedom are utilized. 

Figures 2.9 and 2.10 illustrates the average IR-EE versus the average transmit power, 
and the average achievable rate versus the average transmit power, respectively. The 
curves are obtained by solving Problem 2.4 for co 2 = 0 and 0 < < 1, j = 1,3, where 

the value of coj is uniformly varied with a step size of 0.01 such that coj = 1. Without 
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Figure 2.9.: Average IR-EE versus Eigure 2.10.: Average achievable rate versus 

average transmit power. average transmit power. 


the objective of EH-EE maximization, the power allocation policy in the proposed EE 
algorithm is designed for IR-EE maximization and transmit power minimization. For the 
proposed EE algorithm, we can observe from Figures 2.9 and 2.10 that for small transmit 
power, both IR-EE and the achievable rate grow monotonically as the transmit power 
ascends from zero. Figure 2.9 shows that IR-EE approaches to its maximum point at a 
very small transmit power. Eigure 2.10 shows the corresponding achievable rate which 
remains at a low level due to the small transmit power. In contrast, for the baseline 
scheme, we can see a bell-shaped trend of IR-EE in Figure 2.9 and a monotonically 
ascending trend of the achievable rate in Figure 2.10 with the increasing transmit power. 
In the small transmit power regime, IR-EE and data rate behave similarly as in the EE 
algorithm. However, in the high transmit power regime, the logarithmical growth of 
data rate is slower than the linear increment of the transmit power, which leads to 
energy-inefficient, i.e., IR-EE declines. Besides, when Nj = 8, the system performance 
shows a reduction on IR-EE and an growth on rate compared to the case of 4 transmitting 
antennas. This is because the achievable rate is improved by exploiting extra degrees 
of freedom offered by more transmitting antennas. However, this improvement cannot 
compensate the increment of antenna power consumption. Thus, a lower IR-EE is 
resulted. 

Figures 2.11 and 2.12 depict the average EH-EE versus the average transmit power, 
and average harvested energy versus the average transmit power, respectively. Similarly, 
the curves are obtained by solving the MOOP for = 0 and 0 < coj < l,j = 2,3, where 
the value of cOj is uniformly varied with a step size of 0.01 such that cOj = 1. With no 
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Figure 2.11.: Average EH-EE versus 
average transmit power. 


Eigure 2.12.: Average harvested energy 
versus average transmit 
power. 


concerns on the objective of IR-EE maximization, the power allocation policy is designed 
for EH-EE maximization and transmit power minimization. It can be observed that both 
EH-EE and harvested energy are growing with the increasing transmit power. Especially, 
the curves of the proposed EE algorithm and the baseline scheme overlap, which means 
that the maximal EH-EE and the maximal harvest energy are simultaneously obtained 
by the same amount of transmit power. This is thanks to the linear relationship between 
the harvested power and the transmit power. In terms of the comparison for different 
number of antennas, a better performance on both EH-EE and harvested energy is 
showed for 8 transmitting antennas. In particular, for Nj = 8, we see that EH-EE 
increases faster with the increasing transmit power than the case of = 4. This implies 
a more efficient and effective power transfer is achieved by using more transmitting 


antennas. 
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Chapter 3. 

Power-Efficient SWIPT in Secure 
Communication Systems 

In this chapter, we study resource allocation algorithm design for power-efficient SWIPT 
in secure communication systems with power splitting receivers. In a multiuser system, 
the transmitter supports SWIPT to the desired receivers and power transfer to roaming 
receivers. In particular, the roaming receivers are potential eavesdroppers, thus, artificial 
noise is applied to facilitate secure communication. Under the consideration of system 
QoS, an optimization problem is formulated for minimizing the total transmit power by 
jointly optimizing the beamforming vectors, power splitting ratios at the desired receivers, 
and the power of artificial noise. We proposed a power-efficient resource allocation 
algorithm which enables the dual use of artificial noise for supplying power transfer and 
guaranteeing communication security. The non-convex problem is transformed into SDP 
and solved by SDP relaxation. The global optimum can be achieved. Simulation results 
illustrate a significant power saving via the proposed algorithm. 


3.1. System Model 

We consider a multiuser communication system with SWIPT in downlink scenario. The 
system consists of one transmitter and two types of receivers, namely, desired receivers 
and roaming receivers, cf. Figure 3.1. The transmitter equipped with Nj antennas serves 
K desired receivers and M roaming receivers. The desired receivers are single antenna 
devices with low computational capability. They exploit the received RF signal for both 
information decoding and energy harvesting. On the other hand, the roaming receivers 
are wireless terminals with antennas (Nj > . In particular, they could belong 

to other communication systems and search for additional RF energy supply. Suppose 
they temporally connect to the transmitter for energy harvesting from signals radiated 
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from the transmitter^. However, it is possible that the roaming receivers eavesdrop 
the information signals deliberately. As a result, the roaming receivers are potential 
eavesdroppers which should be taken into account for secure communication. 


a 



Figure 3.1.: Multiuser SWIPT system with power splitting receivers. 


We assume time-slot-based transmission. In each scheduling time slot, K independent 
precoded signal streams are transmitted simultaneously to K desired receivers. Specifi¬ 
cally, a dedicated beamforming vector, e jg assigned to each desired receiver 

to facilitate information transmission. At the same time, the messages intended for the 
desired receivers may be overheard by roaming receivers since all receivers are in the 
range of service coverage. In order to provide secure communication, artificial noise 
is generated and transmitted concurrently to interfere the reception of the roaming 
receivers. As a result, the transmit signal x e jg composed of the K desired 

information signals and the artificial noise, which is given as 

K 

X = ^WfcSfc-Fv, (3.1.1) 

k=i 

where Sj. e C is the signal to the kth desired receiver, k = Without loss 

of generality, we assume = l,Vk. v e c^tx! jg artificial noise vector 

generated by the transmitter to degrade the quality of the received signal of the potential 


possible scenario of the considered system model is a cognitive radio setup. Specifically, the roaming 
receivers may be primary receivers which harvest energy from a secondary transmitter for extending 
the lifetime of the primary network. 
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eavesdroppers. In particular, we model the artificial noise vector as v ~ ‘<^^(0, V) with 
zero mean and covariance matrix V = £{vv^}. 

We focus on a frequency flat fading channel and a TDD system. For the transmitter, 
perfect CSI of all receivers can be obtained by channel reciprocity and handshaking 
signals. The received signals at the desired receivers and the roaming receivers are given 


by 



(3.1.2) 

(3.1.3) 


where h;;. e c^txi denotes the channel vector between the transmitter and desired 
receiver k. The channel matrix between the transmitter and roaming receiver m is 
denoted by e ([;;NTXiVR_ capture the joint effects of multipath fading and 

path loss, ~ and ~ are additive white Gaussian 

noises (AWGN) caused by the thermal noises in the antennas of the desired receivers 
and the roaming receivers, respectively. 

We assume a power splitting structure [6] is adopted in both desired receivers and 
roaming receivers. In particular, at the receiver RF front-end, the received signal at 
desired receiver k is divided into two power streams where 100 x p;.% are used for 
decoding information and the remaining 100 x (1 — pj^)% are used for harvesting 
energy, cf. Figure 3.1. Here, 0 < < 1 is the power splitting ratio of desired receiver 

k. Similarly, power splitting is also performed at the roaming receivers for energy 
harvesting and information decoding. We assume that all receivers have enough energy 
for information decoding at the current time instant which is independent of the amount 
of harvested energy. The harvested energy is stored in battery and used to support 
normal operation of the receiver in the future. Since a portion of received power is 
reserved for energy harvesting, the equivalent received signal model for information 
decoding at desired receiver k can be written as 



(3.1.4) 


where is AWGN with zero mean and variance cr^ caused by signal processing, cf. 
Figure 3.1. We assume that the signal processing noise variances are the same for all 
receivers in this chapter. 
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3.2. Problem Formulation 

In this section, we first introduce the adopted QoS metrics for SWIPT system design. 
Then, the power-efficient resource allocation algorithm design is formulated as a non- 
convex optimization problem. 

The achievable rate (bit/s/Hz) between the transmitter and desired receiver k is given 

by 


Rk = log2(l + rfc), where 

^ ^ _ PfcIhfcWfcp _ 

Pk ( S + Tr(hfch"V) + 0-2^^) + 0-2 

is the SINK at desired receiver k. 

On the other hand, to guarantee communication security, the roaming receivers are 
treated as potential eavesdroppers who attempt to decode the messages transmitted to 
desired receivers. Thereby, we consider the worst case scenario that roaming receivers 
are high computational capable eavesdroppers. We assume that roaming receiver m 
is able to perform successive interference cancellation (SIC) to remove all multiuser 
interference before decoding the signal of receiver k. Therefore, the data rate between 
the transmitter and roaming receiver m for decoding the signal of desired receiver k can 
be represented as 

= log2det(^Ijv^ + A“VJ^G"wfcW^G,„), where (3.2.3) 

= G"VG,„ + (3.2.4) 


(3.2.1) 

(3.2.2) 


0 < < 1 is the power splitting ratio of roaming receiver m. is the interference- 

plus-noise covariance matrix for roaming receiver m. In practice, the roaming receiver 
can be malicious and devote all the received energy for information decoding. Thus, the 
data rate in (3.2.3) is bounded above by 


det + o-sI]vJ 


(3.2.5) 


which is obtained by setting = 1 in (3.2.3). 

Consequently, considering this worst case scenario, the maximum achievable secrecy 
rate of desired receiver k is given by 


R 


sec 

k 



max 



+ 


(3.2.6) 
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Remark 1. We note that the results of this work are also applicable to the case of roaming 
receivers (potential eavesdroppers) employing single user detectors by modifying the term 
in (3.2.4) accordingly. 

In terms of power transfer, both the information signal and the artificial noise serve 
as energy source for the receivers due to wireless broadcasting property. The total 
harvested energy of desired receiver k is given by 

K 

Ek = + (3.2.7) 

J=1 

where 0 < 17 < 1 is the energy conversion efficiency, indicating the efficiency of con¬ 
verting the received RF energy to electrical energy for storage. We assume that it is a 
constant and is identical for all receivers. 

Similarly, the total amount of energy harvested by roaming receiver m is given by 

K 

El = >) J1 - pD (XiTr(G>*<G J + Tr(G„G“V) + . C3.2.8) 

k=l 


The system design objective is to minimize the total transmit power while providing 
system QoS on secure communication and power transfer. The power efficient resource 
allocation algorithm design is formulated as an optimization problem which is given by 


Problem 3.1. 



K 


minimize 

VeH~T,Wfc,pfc 

El 

k=l 

w,|p + Tr(V) 

subject to 

Cl : 

r, > Vfc, 


C2 : 

Reav < T^max 
m,k — m,k ’ 


C3 : 

Ffc > PP\ V/c, 



K 


C4: 

rj(2Tr(G>, 



fc=i 


C5 : 

0 < pfc < 1, V/c, 


C6 : 

VE 0. 


Constraint Cl indicates that SINR at desired receiver k is required to be larger than a 
given threshold, > 0. Since any desired receiver could be chosen as an eavesdropping 
target of roaming receiver m, the upper limit is imposed in C2 to restrict the 
eavesdropping rate of roaming receiver m when it attempts to decode the message of 
desired receiver k. Notice that in practice we are interested in the case of R^. > 
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for secure communication, which means > Ri, — max = loenfl + — 

rnax^ ^ particular, the parameters and can be selected to 

provide flexibility in designing power-efficient resource allocation algorithms for different 
applications. Constants and P^'^^ in constraints C3 and C4 specify the minimum 
required harvested energy at desired receiver k and roaming receiver m, respectively The 
physical meaning of constraint C4 is that the transmitter only guarantees the minimum 
required harvested power at roaming receiver m when it does not intend to eavesdrop, 
i.e., = 0. Constraint C5 implies the physical constraint for the power splitter. In 

addition, we assume that the power splitter is a passive device which does not consume 
or gain any power when splitting the received signal. Constraint C6 and V e ensure 
that the covariance matrix V is a positive semi-definite Hermitian matrix. 


3.3. Power-Efficient Resource Allocation Algorithm 
Design 

It can be observed that Problem 3.1 is non-convex due to constraints Cl and C2. To 
overcome the non-convexity of C2, we recast Problem 3.1 as SDR We first replace Wj^w^ 
by = wpw^ and rewrite C2 as 

C2 : det (ijv^ + Q“^G"WfcG^) < fc, Vm, k, where (3.3.1) 

Q. = G«VG^ + (a-f^, + a-^^)I^^^0. 

is an auxiliary constant that k 'ip,ri,k > 1 for ^ Then, we 

introduce the following proposition to simplify the considered problem. 

Proposition 3.1. For R^^ > 0, Vm, fc, the following implication on constraint C2 holds: 

C2^^: G«WfcG^^('0^,,-l)Q,, Vm,fc. 

C2 is a linear matrix inequality (LMI) constraint. Specifically, constraints C2 and C2 
are equivalent if Rank(W;^) = 1, Vfc. 

Proof. Please refer to Appendix B.3. □ 

Now, we apply Proposition 3.1 to Problem 3.1 by replacing constraint C2 with C2. 
Then, the reformulated problem can be written as 
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Problem 3.2. 


K 


minimize 



subject to Cl : — 

r I- 


k 



C2 : 




K jjreql 


C3 : 


K 


C4: Tr(G«(V + 2w,)G„)> 


Tr(G«(V + 2w,)G„)>^ 

i-=i h 



C5 : 0 < pfc < 1, Vfc, 


C6 : V ^ 0, 


C7 : 0, Vfc, 


C8 : Rank(Wj = 1, Vfc. 


Constraints C7, C8, and e vfc, are imposed to guarantee that holds 

for the optimal solution. In general, replacing constraint C2 by C2 leads to a larger 
feasible solution set for the problem, cf Proposition 3.1. However, Problem 3.1 and 3.2 
are equivalent for Rank(Wjp) = 1, Vk. Thus, we focus on Problem 3.2 in the following. 

Although the new constraint C2 is an affine function with respect to the optimization 
variables, it can be verified that Problem 3.2 is still non-convex due to the combinatorial 
rank constraint in C8. In order to achieve an efficient resource allocation algorithm 
design, we adopt SDP relaxation approach. In particular, we relax constraint C8 by 
removing it from the problem formulation, such that Problem 3.2 becomes a convex 
problem. The SDP relaxed problem is given by 

Problem 3.3. 


K 



(3.3.2) 


We note that Problem 3.3 can be solved efficiently by numerical solvers such as CVX 
[38]. Notably, if the optimal solution of Problem 3.3 admits a rank-one matrix, 
then Problem 3.1, 3.2, and 3.3 share the same optimal solution and the same optimal 
objective value, i.e., the global optimum is achieved. 

Now, we introduce the following theorem to reveal the tightness of the SDP relaxation. 
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Theorem 2. Suppose the optimal solution of Problem 3.3 is denoted by {W^,V*,p^}^ 
> 0^ and >0. If 3k: Rank(wp > 1, then we can construct another solution for 
Problem 3.3, denoted as which not only achieves the same objective value as 

{W^,V*,p^}^ but admits a rank-one matrix, i.e., Rank(Wj^) = 1, Vk. 

Proof. Please refer to Appendix B.4 for a proof of Theorem 2 and a method for construct¬ 
ing {Wj^,V, Pj^} with Rank(W;^) = 1, Vk. □ 

In other words, by applying Theorem 2 and Proposition 3.1, the global optimal 
solution of the original problem is obtained. 


3.4. Results 

In this section, we demonstrate the system performance of the proposed power efficient 
resource allocation design by simulation results. In particular, we solve Problem 3.1 for 
different channel realizations and show the corresponding average system performance. 


Table 3.1.: Simulation Parameters 


Carrier center frequency 

470 MHz 

Number of desired receiver 

K = 3 

Number of roaming receiver 

M = 2 

Number of receiving antenna 

Nr = 2 

Antennas gain 

10 dBi 

Antenna noise power 

aL = -124dBm 

Signal processing noise power 

CT^ = —23 dBm 

S 

Rician factor 

3 dB 

Reference distance 

2 meters 

Maximum service distance 

50 meters 

Minimum required SINR 

preq 

Maximum data rate of roaming receivers 

^ bit/s/Hz 

Minimum required harvested power 

preqi ^ req2 ^ q 
k m 

Energy conversion efficiency 

LD 

d 

11 


The simulation parameters are summarized in Table 3.1. We adopt the TGn path 
loss model [39]. In particular, we assume a carrier center frequency of 470 MHz. 
Assume 3 desired receivers and 2 roaming receivers (potential eavesdroppers), which 
are uniformly distributed in the range between a reference distance of 2 meters and 
a maximum distance of 50 meters. Each roaming receiver is equipped with = 2 
antennas. The multipath fading coefficients are modeled as independent and identically 
distributed Rician fading with Rician factor 3 dB. We set the minimum required SINR 
of all desired receivers identical, i.e., Vk, the maximum tolerable rate of 
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Figure 3.2.: Average total transmit power versus the minimum required SINK of desired 
receivers for Ai-j. = 5 and N-^ — 8. The double-sided arrows indicate the 
power gains achieved by the optimal scheme compared to the baseline 
schemes. 


each roaming receiver is = 1 bit/s/Hz, Vm, k, and the minimum required harvested 
power for all receivers is = 0 dBm. The energy conversion efficiency 17 is 

0.5. Antenna gain is 10 dBi, and the antenna noise power is = —124 dBm at a 
temperature of 290 Kelvin. We assume a 8-bit uniform quantizer is employed in the 
analog-to-digital converter at the analog front-end of each receiver, which result in a 
signal processing noise of cr^ = —23 dBm. 

In Figure 3.2, we show the average total transmit power versus the minimum required 
SINR for different numbers of transmitting antenna and different resource allocation 
schemes. It can be observed that the total transmit power increases monotonically with 
an increasing SINR. The reason behind this is twofold. First, a higher transmit power 
of information signals is required to satisfy the more stringent requirement on SINR. 
Second, a higher amount of power also has to be allocated to the artificial noise to 
neutralize the increasing information leakage potential due to a stronger information 
signal, cf. Figure 3.3. On the other hand, we see that a significant power saving can be 
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Figure 3.3.: Average transmit power allocation for desired signal and artificial noise 
versus the minimum required SINK of desired receivers for N-^ = 8. 

achieved by the proposed optimal scheme when the number of antenna increase from 
Nj = 5 to Nj' = 8. This is because more transmitting antennas provide extra degrees of 
freedom, which enables a more power-efficient resource allocation. 

For comparison, we also present the performance of two simple suboptimal baseline 
schemes. For baseline scheme 1, zero-forcing beamforming is applied for information 
signal such that the desired receivers do not experience any multiuser interference. In 
particular, we calculate the eigenvalue decomposition of for desired 

receiver k where H.;;. = [h^ ... hj^+i ... is an Nj x Nj unitary matrix, and 

is a diagonal matrix with ascending eigenvalues of as main diagonal elements. 

Then, we select W^. = qsu\'^suhi,'^suhp where > 0 is a new scalar optimization 
variable and is the first column vector^ of U;;. such that = 0. In other 

words, the directions of the beamforming matrices are fixed for all desired receives. 
Then, we minimize the total transmit power by optimizing qsu\,^’ subject to the 

^In general, different column vectors with respect to the null space of can be used as zero¬ 

forcing beamforming vector. For algorithm computational simplicity, we select the first column vector 
corresponding to the minimum eigenvalue of matrix as zero-forcing beamforming vector. 
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Figure 3.4.: Average secrecy rate per desired receiver versus the minimum required SINR 
of desired receivers. 


same constraints as in 3.3. We note that the zero-forcing beamforming matrix admits 
a rank-one structure. As for baseline scheme 2, it shares the same resource allocation 
policy as baseline scheme 1 except that we set pj. = 0.5, Vk. It can be observed in 
Figure 3.2 that the optimal scheme achieves significant power savings over the two 
baseline schemes. Notably, the performance gain of the optimal scheme over the two 
baseline schemes is further enlarged by more transmitting antennas. Since the optimal 
scheme can fully utilize the degrees of freedom for resource allocation. In contrast, 
although multiuser interference is eliminated in the two baseline schemes, the degrees 
of freedom for resource allocation in the baseline schemes are limited thus resulting in a 
higher transmit power. Furthermore, the performance gap between baseline scheme 1 
and baseline scheme 2 reveals the performance gain on the optimization of the power 
splitting ratio. 

Figure 3.3 depicts the average transmit power allocated to the information signal and 
the artificial noise for = 8. As shown, the power allocations for both information 
signal and artificial noise grow rapidly with the increment of minimum SINR requirement. 
Besides, both the optimal scheme and the two baseline schemes indicate that a large 
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Figure 3.5.: Average total harvested power versus the minimum required SINK of desired 
receivers. 


portion of the total transmit power is allocated to the artificial noise. These results 
suggest that artificial noise generation is crucial for supporting communication security 
and WPT. 

Figure 3.4 plots the average secrecy rate per desired receiver with respect to the 
minimum required SINR of desired receivers for different numbers of transmitting 
antenna and different resource allocation schemes. The average achievable secrecy rate 
rises with the minimum required SINR since the eavesdropping rate of roaming receivers 
(potential eavesdroppers) is limited toR““ = 1 bit/s/Hz. Besides, all considered schemes 
are able to guarantee the QoS requirement on communication security (constraints Cl 
and C2) and achieve the same secrecy rate. However, the two baseline schemes achieve 
the same secrecy rate as the optimal scheme at the expense of a significantly higher 
transmit power, cf. Figure 3.2. 

In Figure 3.5, we demonstrate the average total harvested power versus the minimum 
required SINR of desired receivers for different resource allocation schemes and different 
numbers of transmitting antenna. The average total harvested power is computed by 
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assuming the roaming receivers (potential eavesdroppers) do not eavesdrop. It is shown 
that the average total harvested power of all resource allocation schemes increase with 
As more RF energy is available due to stronger information signals, cf. Figure 
3.2. Besides, more energy is harvested in the two baseline schemes compared to the 
optimal scheme. The superior energy harvesting performance of the baseline schemes 
comes at the expense of a significant large transmit power. On the other hand, it can 
be observed that the average total harvested power in the system decreases with an 
increasing number of transmitting antenna. This is because when extra degrees of 
freedom are offered by more transmitting antennas, system resource allocation becomes 
more efficient. In other words, the information leakage can be efficiently reduced since 
artificial noise jamming can be more accurately performed. As a results, a lower transmit 
power is required to fulfill the considered QoS requirements and less power is harvested. 



38 



39 


Chapter 4. 
Conclusion 


In this thesis, we focused on energy-efficient resource allocation algorithm design for 
SWIFT and considered the secure communication in SWIFT systems. 

In Chapter 2, we investigated the energy efficiencies for information transmission and 
energy harvesting in a basic SWIFT system with separated receivers, i.e., one information 
receiver and one energy harvester. The trade-off between three desirable but conflicting 
system design objectives, namely, IR-EE maximization, EH-EE maximization, and total 
transmit power minimization, was studied by formulating a MOOE The non-convex 
problem was transformed by Charnes-Cooper transformation method and solved by SDF 
relaxation approach. The complete Fareto optimal set of the MOOF was achieved by 
the proposed energy-efficient resource allocation algorithm. In particular, the algorithm 
provided flexibility in power allocation when balancing between multiple objectives. 
Simulation results showed the trade-off region of the considered objectives and revealed 
the system performance gain on energy efficiency compared to the baseline scheme. 

In Chapter 3, we studied the power-efficient resource allocation algorithm design 
for secure communication in a SWIFT system with power splitting receivers. Under 
the eavesdropping potential of multi-antenna energy harvesters, we adopted artificial 
noise generation to guarantee secure information delivery to legitimate receivers. The 
algorithm design was formulated as a non-convex optimization problem taking into 
account the QoS requirements on communication security and efficient power transfer. 
We applied SDF relaxation approach to obtain the optimal solution. Especially, the 
proposed power-efficient algorithm achieved the dual use of artificial noise on ensuring 
communication security and facilitating EH. Simulation results confirmed the remarkable 
performance of the proposed optimal scheme on energy saving and communication 
security. 

In the future work, we are interested in energy-efficient resource allocation in orthog¬ 
onal frequency multiple access (OFDMA) systems with secure SWIFT. 



40 



Bibliography 


[1] Powercast Coporation, “RF Energy Harvesting and Wireless Power for Low-Power 
Applications,” 2011. [Online], Available: http://www.mouser.com/pdfdocs/ 
Powercast- Overview-2011-01 -25 .pdf 

[2] I. Krikidis, S. Timotheou, S. Nikolaou, G. Zheng, D. W. K. Ng, and R. Schober, 
“Simultaneous Wireless Information and Power Transfer in Modern Communication 
Systems,”/£££ Commun. Mag., vol. 52, no. 11, pp. 104-110, Nov. 2014. 

[3] Z. Ding, C. Zhong, D. W. K. Ng, M. Peng, H. A. Suraweera, R. Schober, and 
H. V Poor, “Application of smart Antenna Technologies in Simultaneous Wireless 
Information and Power Transfer,” 2015, to appear in the IEEE Commun. Mag. 

[4] L. Varshney “Transporting Information and Energy Simultaneously,” in Proc. IEEE 
Intern. Sympos. on Inf. Theory, Jul. 2008, pp. 1612 -1616. 

[5] P Grover and A. Sahai, “Shannon Meets Tesla: Wireless Information and Power 
Transfer,” in Proc. IEEE Intern. Sympos. on Inf Theory, Jun. 2010, pp. 2363 -2367. 

[6] X. Zhou, R. Zhang, and C. K. Ho, “Wireless Information and Power Transfer: 
Architecture Design and Rate-Energy Tradeoff,” IEEE Trans. Commun., vol. 61, pp. 
4754-4767, Nov. 2013. 

[7] R. Zhang and C. K. Ho, “MIMO Broadcasting for Simultaneous Wireless Information 
and Power Transfer,” in Proc. IEEE Global Telecommun. Conf, Dec. 2011, pp. 1 -5. 

[8] J. Xu, L. Liu, and R. Zhang, “Mutiuser MISO Beamforming for Simultaneous 
Wireless Information and Power Transfer,” IEEE Trans. Signal Process., vol. 62, pp. 
4798 - 4810, Jul. 2014. 

[9] Q. Shi, L. Liu, W. Xu, and R. Zhang, “Joint Transmit Beamforming and Receive 
Power Splitting for MISO SWIPT Systems,” submittedfor possible journal publication, 
Sep. 2013. 



42 


Bibliography 


[10] J. Xu, L. Liu, and R. Zhang, “Multiuser Beamforming for Simultaneous Wireless 
Information and Power Transfer,” Proc. IEEE Intern. Conf. on Acoustics, Speech and 
Signal Process., pp. 4754-4758, May 2013. 

[11] D. W. K. Ng, E. S. Lo, and R. Schober, “Energy-Efficient Resource Allocation in 
Multiuser OFDM Systems with Wireless Information and Power Transfer,” in Proc. 
IEEE Wireless Commun. and Networking Conf., 2013. 

[12] D. W. K. Ng and R. Schober, “Spectral Efficient Optimization in OFDM Systems 
With Wireless Information and Power Transfer,” in Proc. Europ. Signal Process. 
Conf, Sep. 2013, pp. 1-5. 

[13] D. W. K. Ng, E. S. Lo, and R. Schober, “Wireless Information and Power Transfer: 
Energy Efficiency Optimization in OFDMA Systems,” IEEE Trans. Wireless Commun., 
vol. 12, pp. 6352-6370, Dec. 2013. 

[14] J. Park and B. Clerckx, “Transmission strategies for joint wireless information and 
energy transfer in a two-user MIMO interference channel,” in Proc. IEEE Intern. 
Commun. Conf, 2013, pp. 591 -595. 

[15] R. Morsi, D. Michalopoulos, and R. Schober, “Multi-User Scheduling Schemes 
for Simultaneous Wireless Information and Power Transfer,” in Proc. IEEE Intern. 
Commun. Conf, Jun., pp. 4994—4999. 

[16] D. W. K. Ng and R. Schober, “Resource Allocation for Coordinated Multipoint 
Networks With Wireless Information and Power Transfer,” in Proc. IEEE Global 
Telecommun. Conf, Dec. 2014, pp. 4281-4287. 

[17] M. Ch 5 monova, R. Morsi, D. W. K. Ng, and R. Schober, “Optimal Multiuser 
Scheduling Schemes for Simultaneous Wireless Information and Power Transfer,” 
Feb. 2015. [Online]. Available: http://arxiv.org/abs/1502.02179 

[18] Q. Wu, M. Tao, D. W. K. Ng, W. Chen, and R. Schober, “Energy-Efficient Trans¬ 
mission for Wireless Powered Multiuser Communication Networks,” in Proc. IEEE 
Intern. Commun. Conf, Jun. 2015. 

[19] A. A. Nash, X. Zhou, S. Durrani, and R. A. Kennedy, “Wireless Energy Harvesting 
and Information Relaying: Adaptive Time-Switching Protocols and Throughput 
Analysis,” submitted for possible journal publication, Oct. 2013. [Online]. Available: 
http://arxiv.org/abs/1310.7648 




Bibliography 


43 


[20] G. Zheng, Z. Ho, E. A. Jorswieck, and B. Ottersten, “Information and Energy 
Cooperation in Cognitive Radio Networks,” IEEE Trans. Signal Process., vol. 62, pp. 
2290 - 2303, Mar. 2014. 

[21] A. D. Wyner, “The Wire-Tap Channel,” Tech. Rep., Oct. 1975. 

[22] S. Goel and R. Negi, “Guaranteeing Secrecy using Artificial Noise,” IEEE Trans. 
Wireless Commun., vol. 7, pp. 2180 - 2189, Jun. 2008. 

[23] Q. Li and W. K. Ma, “Spatically Selective Artificial-Noise Aided Transmit Optimiza¬ 
tion for MISO Multi-Eves Secrecy Rate Maximization,” IEEE Trans. Signal Process., 
vol. 61, pp. 2704-2717, May 2013. 

[24] D. W. K. Ng, E. S. Lo, and R. Schober, “Energy-Efficient Resource Allocation for 
Secure OFDMA Systems,” IEEE Trans. Veh. TechnoL, vol. 61, pp. 2572 - 2585, May 
2012 . 

[25] -, “Secure Resource Allocation and Scheduling for OEDMA Decode-and- 

Eorward Relay Networks,” IEEE Trans. Wireless Commun., vol. 10, pp. 3528 - 
3540, Aug. 2011. 

[26] D. W. K. Ng, M. Shaqfeh, R. Schober, and H. Alnuweiri, “Robust Layered 
Transmission in Secure MISO Multiuser Unicast Cognitive Radio Systems,” 
submitted for possible journal publication, Jun. 2014. [Online]. Available: 
http: / / arxiv. org/abs/1406.6542 

[27] L. Liu, R. Zhang, and K. C. Chua, “Secrecy Wireless Information and Power Transfer 
with MISO Beamforming,” IEEE Trans. Signal Process., vol. 62, pp. 1850 - 1863, 
Jan. 2014. 

[28] D. W. K. Ng, E. S. Lo, and R. Schober, “Robust Beamforming for Secure 
Communication in Systems with Wireless Information and Power Transfer,” 
accepted for publication, IEEE Trans. Wireless Commun., Mar. 2014. [Online]. 
Available: http://arxiv.org/abs/1311.2507 

[29] D. W. K. Ng., L. Xiang, and R. Schober, “Multi-Objective Beamforming for Secure 
Communication in Systems with Wireless Information and Power Transfer,” in Proc. 
IEEE Personal, Indoor and Mobile Radio Commun. Sympos., Sep. 2013, pp. 8-12. 

[30] D. W. K. Ng and R. Schober, “Resource Allocation for Secure Communication 
in Systems with Wireless Information and Power Transfer,” in Proc. IEEE Global 
Telecommun. Conf, Dec. 2013. 




44 


Bibliography 


[31] S. Leng, D. W. K. Ng, and R. Schober, “Power Efficient and Secure Multiuser 
Communication Systems with Wireless Information and Power Transfer,” in Proc. 
IEEE Intern. Commun. Conf., Jun. 2014, pp. 800-806. 

[32] D. W. K. Ng, R. Schober, and H. Alnuweiri, “Secure Layered Transmission in 
Multicast Systems With Wireless Information and Power Transfer,” in Proc. IEEE 
Intern. Commun. Conf., Jun. 2014, pp. 5389-5395. 

[33] D. W. K. Ng and R. Schober, “Secure and Green SWIPT in Distributed Antenna 
Networks with Limited Backhaul Capacity,” submitted for possible journal 
publication, Oct. 2014. [Online]. Available: http://arxiv.org/abs/1410.3065 

[34] -, “Max-Min Fair Wireless Energy Transfer for Secure Multiuser Communication 

Systems,” in Proc. IEEE Inf Theory Workshop, Nov. 2014, pp. 326-330. 

[35] D. W. K. Ng, E. S. Lo, and R. Schober, “Multi-Objective Resource Allocation for 
Secure Communication in Cognitive Radio Networks with Wireless Information 
and Power Transfer,” submitted for possible journal publication. Mar. 2014. 
[Online]. Available: http://arxiv.org/abs/1403.0054 

[36] R. T. Marler and J. S. Arora, “Survey of Multi-objective Optimization Methods for 
Engineering,” Structural and Multidisciplinary Optimization, vol. 26, pp. 369-395, 
Apr. 2004. 

[37] S. Boyd and L. Vandenberghe, Convex Optimization. Cambridge University Press, 
2004. 

[38] M. Grant and S. Boyd, “CVX: Matlab software for disciplined convex programming, 
version 2.0 beta,” [Online] https://cvxr.com/cvx, Sep. 2012. 

[39] IEEE P802.ll Wireless LANs, “TGn Channel Models", IEEE 802.11-03/940r4, Tech. 
Rep., May 2004. 

[40] A. Charnes and W. W. Cooper, “Programming with Linear Fractional Functions,” 
Naval Res. Logist. Quart, vol. 9, pp. 181-186, Apr. 1962. 




45 


Appendix A. 

Mathematical Preliminaries 


A.l. Optimization Problem 


A typical optimization problem consists of an objective function, constraints, and op¬ 
timization variables. In general, the objective function is a function of optimization 
variables that evaluates one aspect of system performance. The set of constraints map the 
considered QoS requirements and physical limitations of the system. A single-objective 
optimization problem is given standardly by 

minimize /o(x) (A. 1.1) 

X 

subject to g;(x) < 0, 1 = 1,2, ...,L, 

h„(x) = 0, n = l,2,...,N, 

where /o(0 is a system objective function and x is a vector of optimization variables. L 
and N are the numbers of inequality constraints and equality constraints, respectively. 
g;(x) is the 1-th inequality constraint and li„(x) is the n-th equality constraint. 

However, in practice multiple desirable system design objectives arise naturally in 
resource allocation problems. As the objective functions could be conflicting with each 
other, non-trivial trade-off occurs in this case, where the solution of single-objective 
resource allocation may not result in satisfactory system performance. Therefore, MOO 
is applied to address this type of resource allocation problem. A standard form of a 
MOOP can be posed as follows [36] : 

minimize F(x) = [Fi(x),F 2 (x),.. .,Fjf(x)]^ (A.l.2) 

X 

subject to g;(x) < 0, 1 = 1,2, ...,L, 

h„(x) = 0, n = l,2,...,N, 
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where K is the number of objective functions and F^.(x),'ik e fC}, is the fc-th 

objective function. 

In contrast to single-objective optimization, a solution to a MOOP is more of an 
abstract concept than a fixed point. In general, there is no single global solution which 
optimizes all the objective functions simultaneously. T 5 ^ically, it is often necessary to 
determine a set of points that fit a predetermined definition for a optimum, which is 
Pareto optimality. Pareto optimality of MOOP is defined as 

Definition A.I. Pareto Optimal: A point, x* G is Pareto Optimal if and only if there does 
not exist another point, xG such that F(x) :< F(x*) and Ffc(x) < Fj^(x*), k = 
for at least one function. 

Evidently, any point that is not in the Pareto optimal set is strictly suboptimal. The 
Pareto optimal set is an analogy to global optimality that can achieve in MOO. We note 
that single-objective optimization problems are special case of MOOPs with K = 1. In 
other words, if a resource allocation algorithm can solve the MOOI^ then it can be used 
to solve the corresponding single-objective optimization problem. 
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Calculations 


B.l. Proof of Proposition 2.1 


The proof is based on the Charnes-Cooper transformation [40] and follows a similar 
approach in [29]. By substitute the new optimization variables in (2.3.1) into the 
original problem 2.1, we can obtain an equivalent problem representation as 


maximize 

w,,WEe]H~T,e 

subject to 


01og2(l + 


= 


Tr(HW,) 

Oaf 


) 


Tr(Wi+WE) 




a, : 0 > 0. 


(B.1.1) 


Now we show that the above Problem (B.1.1) is equivalent to Problem 2.5. First, it 
can be observed that in Problem (B.1.1) 0 = 0 is impossible. Otherwise, Wj = = 0 

according to Cl and C2, the objective function is invalid. Thus, without loss of generality, 
the constraint 0 > 0 can be replaced by 0 >0. Second, we prove by contradiction that 
C4 in Problem 2.5 is satisfied with equality for the optimal solution. Denote the optimal 
solution of Problem 2.5 as (Wj,Wg, 0*). Then, 

Tr(W*+W*) 

-^^ + PJ = 1. (B.1.2) 

Assume that C4 is fulfilled with strict inequality for the optimal solution, i.e., _|_ 

0*(Ai2'Pant+Fc) < 1- Then, we construct a new feasible solution (Wj, W^,, 9') = (cWj, W^, c0*), 
where c > 1, such that + P^) = 1. It can be verihed that (Wj, Wg, 6') 

achieves a larger objective value in Problem 2.5 than (Wj, W^, 0*). Then, (Wj, W^, 0*) 
cannot be the optimal solution. Contradiction occurs. Thus,C4 must hold with equality. 

The equivalency between Problem (B.1.1) and Problem 2.5 is proved, which means 
Problem 2.5 is equivalent to the original Problem 2.1. Similarly, The equivalency of 
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Problem 2.6, 2.7 and 2.8 to their original problems can be proved by following the above 
approach. 


B.2. Proof of Theorem 1 


Theorem 1 can be proved by analyzing the KKT optimality conditions of the SDP relaxed 
Problem 2.9. First we need the Lagrangian function as the following 


if (Wi, We, 6 , T, a, X, Y, n, 72, Ts, 5) 

T + a(Tr(Wi + We) - - Tr(XWi) - Tr(YWE) -56 


(B.2.1) 


Tr(Wi + WE) , r . ^ Tr(HWi) , 

+ Pi --1" ^(•^T^’ant + ^’c) “ 1) +Tl -—)) 


+ 72 


(1 - ^ Tr(G(Wi + We))) -t]+ 7 . 


COr, 


IR 

? .1 


T 


( Q ^T^ant Pc 


)-t]. 


where a, ^,X, Y, 72 , 73 ,5 are dual variables associated with the corresponding con¬ 
straints, respectively. Since the SDP relaxed Problem 2.9 satisfies Slater’s constraint 
qualification and is convex with respect to the optimization variables, strong duality 
holds. Denote the optimal solution as {Wj, We, 6 *, t*}, and the optimal dual variables as 
{a*, 7 *, y*, y*, 5*}. Then, based on KKT optimality conditions, the gradient of 

Lagrangian function with respect to Wj and We vanish, and the complementary slackness 
condition is satisfied as well. Thus, we have 


6 

awj ? 


X 


yi(Ui0 


d S£ B 72^2V 

^ = (a + -)I-Y-^^G 
5We_ ? 

Tr(XWi) = 0, 


4>*j,(0o-f + Tr(HWi))' 
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-H 


72^2V 

EH 


G = 0 


^EH 


Tr(YWE) = 0. 

P 

X=ia+'-)l- 

P 

Y=(a + -)I- 


72^2V 

^EH 


G- 


EH 


$*^(0o-2 + Tr(HWi)) 


H 


XWi = 0, 
YWe = 0. 


(B.2.2) 

(B.2.3) 

(B.2.4) 

(B.2.5) 


Now, we investigate W; and We in three cases. 



















B.3. Proof of Proposition 3.1 
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Case 1: co^ = C 02 = 0, co^ = 1. (B.2.2) and (B.2.3) become X = Y = (a + |)I, which 
means that X and Y are matrices of full rank. Due to (B.2.4) and (B.2.5), Wj lies in the 
null space of X, and Wg lies in the null space of Y. Thus, Wj = Wg = 0. 

Case 2: coi = 0, CO 2 ^ 0. This means that IR-EE maximization is not considered 
and EH-EE maximization is required. Thus, Y 2 ^ ^ must hold to keep constraint C6 
for j = 2. Then, (B.2.2) and (B.2.3) become X = Y = (a + j)I — Since 

^ ^EH ^EH 

is of rank one, based on the basic property of the rank of matrices, Rank(X) > — 1 

and Rank(Y) > Ai-p — 1. If Rank(X) = Rank(Y) = Alp, then Wj = We = 0. which 
can not occur since EH-EE maximization is required. As a result, we can conclude 
Rank(X) = Rank(Y) = Alp — 1, Rank(Wi) = 1 and Rank(WE) = 1. 

In particular, we can see from the problem formulation that Wj and Wg are equivalent 
optimization variables in this case. Thus, the optimization variables can be redefined as 


Wj = Wj + Wg and W^ = 0. Then, the Lagrangian multiplier X^ with respect to Wj results 
-''T — aforementioned, Rank(X') = A/p — 1 and Rank(Wj) = 1 




inX^ = (a + |)I 
must hold. 

Case 3: coj ^ 0, CO 2 ^ 0. IR-EE maximization and EH-EE maximization are both 
considered, thus constraint C 6 for j = 1,2 is active, i.e., Ti 7 ^ 0 and Y 2 7 ^ 0. Since the 
Lagrangian multiplier Y ^ 0, we have a + j according to equation (B.2.3), 

where gj is the largest eigenvalue of G [ 8 ] . If the equality holds, i.e., a + | = -^-^^gi, 

^ "^EH 

we consider X in (B.2.2). Since 0 > 0 is proved in Appendix B.l, to satisfy X ^ 0, = 0 

must hold. Contradiction occurs. This implies a + j > Thus, Y 0, which 

_ ^ "^EH 

means Y is a full-rank matrix. Therefore, Wg = 0. On the other hand, according to 


equation (B.2.2), X = Y — 


yi<^i 6 


-H. Since Rank 


yi<^i 6 


$*^(0crf+Tr(HW,)) 


H 


: 1 , we 
1 must 


<&*^(0a2+Tr(HW,)) 

have Rank(X) > A/p — 1. As aforementioned, Rank(X) = A/p — 1 and Rank(Wi) 
hold. 

Consequently, in all cases, the optimal solution of the relaxed Problem 2.9 satisfies 


Rank(Wj) = 1 and Rank(Wg) < 1. In particular, an optimal solution with Rank(Wj) = 1 
and Wg = 0 can always be constructed. 


B.3. Proof of Proposition 3.1 

(B.3.1) 
(B.3.2) 


We start the proof by expressing constraint C2 as 


(a) 


det(Ijv, + Q;i/2G>,G,„Q;'/^) < 
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where (a) is due to the fact that det(I + AB) = det(I + BA) holds for any matrices A and 
B. Then, we introduce the following lemma which provides a lower bound on the left 
hand side of (B.3.2). 


Lemma 1. For any square matrix A ^ 0, we have det(I + A) > 1 + Tr(A) [23], where the 
equality holds if and only i/Rank(A) < 1. 


Exploiting Lemma 1, the left hand side of (B.3.2) is bounded below by 

det(I^^ + > 1 + Tr(Q;'/2G>,G^Q;'/2). (B.3.3) 

Subsequently, by combining equations (B.3.1), (B.3.2), and (B.3.3), we have the follow¬ 
ing implications: 


(B.3.1) « (B.3.2) (B.3.4a) 

TriQ-J^^Gyj,G„,Q-J/^) < - 1 (B.3.4b) 

- 1 (B.3.4C) 

« G"W,G,„ ^ (t/.^,, - 1)Q^. (B.3.4e) 


Amax(A) denotes the maximum eigenvalue of matrix A. (b) is due to Tr(A) > Ajjiax(A) for 
a positive semidefinite matrix A h 0. We note that equations (B.3.1) and (B.3.4e) are 
equivalent when Rank(W;i.) = 1, Vk. 


B.4. Proof of Theorem 2 


We follow a similar approach as in [27, 28] to prove Theorem 1. The proof is divided 
into two parts. In the first part, we study the solution structure of Problem 3.3. Then 
in the second part, we propose a simple method for constructing an optimal solution 
with rank-one In order to verify the tightness of the adopted SDP relaxation, we 
analyze the Karush-Kuhn-Tucker (KKT) conditions of the SDP relaxed Problem 3.3 by 
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introducing the corresponding Lagrangian and the dual problem. The Lagrangian of 
Problem 3.3 is given by 


(y^k’ Pfcj ^m,k’ ^k’ Pk’ ^m) 

K K 

= 2]Tr(W,) + Tr(V)-Tr(YV)-2Tr(Z,W,) 


(B.4.1) 


k=l 

K 


fc=l 

K 


+ Xi t - psq Tr(hfchfWfc) + X Tr(hfch"w^.) + Tr(hfch"V) + + —a]] 


fc=i 

K 


+ X^fc 

k=l 

M 

^ z 


k 

jreql 




Pk 


k=l 

M preq2 

V 

m 

m=l 
M K 


K 

-WRCTf„-Tr(G„G«(V + X;wJ) 


V 


k=l 




m=lk=l 


where X^j^, Y, and are the dual variable matrices for constraints C2 , C6, and C7, 
respectively, a^, and are the scalar dual variables of constraints Cl, C3, and C4, 
respectively. On the other hand, constraint C5 for pj. is satisfied automatically and the 
optimal pp. will be illustrated in the later part of this proof. Then, the dual problem of 
the SDP relaxed Problem 3.3 is given by 


maximize minimize ^ CWj^,'V, pj^,Zj^,Y,X^,^,a^, 

ak,Pk,'Vm>0 Pfc,Wfc,V£H~T 


(B.4.2) 


Since Problem 3.3 satisfies Slater’s constraint qualification and is jointly convex with 
respect to optimization variables, strong duality holds. Thus solving B.4.2 is equiva¬ 
lent to solving Problem 3.3. We define {W^,V*,p*} and {Z*,Y*,X^ j., v*,^*, a^} as the 
optimal primal solution and the optimal dual solution of Problem 3.3 with Z^,X* 

0, a^, > 0. Now, we focus on those KKT conditions which are useful in the proof: 


where 


Kn = 0, 


a 


Z*=U,-(^* + ^)h,h", 


M 


= k + S + Xi(«; - /5;)hjh* 

;/fc 


m=l 


Pk = 




k^lp 




reql 

k 


(B.4.3) 

(B.4.4) 

(B.4.5) 


(B.4.6) 
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It can be observed from (B.4.6) that constraint C5 is automatically satisfied. Besides, 
“fc’ Pk ^ ^ must holds for > 0 and > 0. On the other hand, because of the 
complementary slackness condition on in (B.4.3), lies in the null space of for 
W* / 0. In other words, the structure of depends on the space spanned by Z^. Thus, 
we focus on the following two cases to reveal the space spanned by Z*^. Without loss of 
generality, we denote = Rank(Uj,). In the hrst case, we investigate the structure of 
when Ufc is a full-rank matrix, i.e., = Afr. By exploiting (B.4.4) and a basic inequality 

for the rank of matrices, we have 

* 

Rank(Zp + Rank((^ + > Rank(Uj 

^ k 

Rank(Zp>NT-l for a*,^*>0. (B.4.7) 


For > 0 and Rank(U^.) = Nj, Rank(wp = 1 and Rank(zp = Nj — 1 must hold 
simultaneously. Next, we consider the case when Rank(Ufc) is rank-deficient, i.e., r^. < N-^. 
Without loss of generality, we define Null(Ufc) = ^ cWtxCWt-u) U^N;;. = 0 

and Rank(Nj,) = Nj — r^. Let s 1 < tj^ < Nj — r^, denote the column 

vector of N^. Then, by exploiting (B.4.4), we have the following equality: 


o^Z* o = —(_ — 

^ k 

Combining Z* ^ 0 and ^ > 0, 

holds in (B.4.8). In other words. 


+ (B.4.8) 

+ ^*)e«h,h«e,^ = 0, Vt, e {1,...,iVT - rJ, 


Z*Nfc = 0 and h^h^N^ = 0 


(B.4.9) 


hold and N;;. lies in the null spaces of and Z'^ simultaneously. Furthermore, 
Rank (Nu 11(Z*)) > Afr — Tj, holds for satisfying Z^N^ = 0. On the other hand, from 
(B.4.7) and Rank(Ufc) = tj^, we obtain 


Rank(Zp > Tfc - 1. (B.4.10) 

Then, by utilizing (B.4.9) and (B.4.10), Rank (Null(zp) is bounded between 

AIt - tfc + 1 > Rank (Null(Zp) > Nj - r^. (B.4.11) 

As a result, either Rank ( Null(zp) = Nj — r^. or Rank ( Null(zp) = Afr — -I-1 holds for 

the optimal solution. Suppose Rank (Null(zp) = Al-j. — r,, and thus Null(zp = Nj^.. Then, 
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we can express as W* = some positive constants >0,yt,^G 

{1,...,N-^ — r^}. Yet, due to (B.4.9), 


Afp-rj. 

Tr(hkhX) = E r.JrCefh^hfe J =o (B.4.12) 

tfc=l 

holds which cannot satisfy constraint Cl for > 0. Thus, Rank ( Null(zp) = Nj—ri^+1 
has to hold for the optimal W^. Besides, there exists one subspace spanned by a unit 
norm vector U;;. e c^txi such that Z*U;. = 0 and N^Uji, = 0. Therefore, the orthonormal 
null space of Z^ can be presented as 

Null(Z*) = {NfcUUfc}. (B.4.13) 

In summary, without loss of generality, we can express the optimal solution of as 

JVp-rfc 

w* = 2] Yt.Qt.Ql + AUfcu" (B.4.14) 

tfe=i 

where /;;. > 0 is some positive scaling constant.In the second part of the proof, for 
Rank(wp > 1, we reconstruct another solution of the Problem 3.3, {W;.,V, p;^}, based 
on (B.4.14) .Let the constructed solution set be given by 

^ Np-rfc 

Wfc = Au,u" = w*-2] Tt.et.e", 

tfc=i 

Wx-n 

V = v*+^ Pk = pl 

tfc=i 

It can be easily verified that {Wji.,V, p,,} not only satisfies the constraints in Problem 3.3, 
but also achieves the same optimal objective value as {W;^, V, p;.} with Rank(Wfc) = 1, Vk. 
The actual values of can be obtained by substituting (B.4.15) and (B.4.16) 

into Problem 3.3 and solving the resulting convex optimization problem for /jj. and 


(B.4.15) 

(B.4.16) 



